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Abstract. 

In this paper we revisit the "eccentric disc instability" , an instability which 
occurs in coherently eccentric discs of stars orbiting massive black holes (MBHs) 
embedded in stellar clusters, which results in stars achieving either very high or 
low eccentricities. The preference for stars to attain higher or lower eccentrici- 
ties depends significantly on the density distribution of the surrounding stellar 
cluster. Here we discuss its mechanism and the implications for the Galactic 
Centre, home to at least one circum-MBH stellar disc. 



1. Introduction 



The Galactic Centre (GC) hosts at least one clockwise rotating disc of O and 
Wolf-Rayet stars, most likely a result of the gravit ational fragmentation of a 
gaseous accretion disc ([Levin &; Beloborodovl 20031 ). These stars are located 
at projected distan c es of 0.05 — 0.5 pc, and have an ave r age age of ~ 6 Mv r 
(jGenzel et al.1 1200.1 : iPaumard et~aXl l200fil : iLu et all l200flt IBartko et all 12003 ). 
Recent hydro dynamical simulations suggest that young stars forming in such 
discs can retain sim i lar in itial orbital eccentricities, which are often non-zero 
(e.g. Bonnell Rice! 120081 ) . The dynamics of stars in such eccentric discs can 
be ver y different from circ ular structures. In particular there exists an insta- 
bility (iMadigan etaD l200fll ) that occurs in coherently eccentric discs which are 
embedded in nuclear stellar clusters such as the one in our GC (see review by 
Schodel in this volume). This instability can propel stars to very high, or low, 
orbital eccentricities within ~ 1 Myr, overcoming the long relaxational time 
scales associated with galactic nuclei. 

The background stellar cluster (or cusp) is an integral part of the instabil- 
ity. In this paper we identify the reasons for this, vary the density profile of 
the cluster to see the effect it has on the instability, and address the possible 
implications of the observed flattened density profil e of the cusp of late- type 
stars in the GC (jBuchholz et al.ll2009l : iDo et alJl2009l : IBartko et alll2009l ). 



2. The Eccentric Disc Instability 

To describe the physics of the instability we work with a simplified template of 
a galactic nucleus. We consider a disc of stars with mass M^ sc orbiting a MBH 
of mass M.. Both are embedded in a power-law stellar cusp of mass M cusp . 
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We make the assumptions that (1) Md; sc <C M cusp , such that the apisidal 
precession of stellar orbits within the disc is driven by the cusp, and (2) initially 
each stellar eccentricity vector (which points to periapse) is aligned and similar 
in magnitude to those of the other stars in the disc. The stellar orbits precess 
with retrograde motion due to the presence of the cusp, i.e., in the direction 
opposite to the orbital rotation of the stars. The precession time scales as 

tprcc- AT/ M \ P(a)f(e) cxa^ 2 f(e), (1) 



where P(a) = lit \J a? j GM* is the period of a star with semi-major axis a, 
N(< a) is the number of stars in the cusp within a, m is the individual mass of 
the stars, and 7 is the power-law index for the space-density profile of the cusp 
p(r) oc r~ 7 . For simplicity let us assume a power-law index 7 = 3/2, such that 
tprec is constant for all values of a. The key element of the instability is that 
/(e) in Equation ([TJ is an increasing function of eccentricity, df /de > 0. 
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Figure 1. This schematic shows the stellar disc from above its orbital plane. 
The stars move counterclockwise on their orbits and precess in the opposite 
direction (dashed arrow). A "test star" is also shown, with a more eccentric 
orbit, and lagging behind in precession. The direction of its velocity vector, 
v, at apoapse is opposite to the coherent gravitational force, F, it feels from 
the rest of the stars in the disc. This causes the torque, r = r x F, acting on 
the star from the disc to decrease its angular momentum, J = r x v. 

The instability works as follows: A star that has a slightly more eccentric 
orbit than the average star in the disc consequently has a larger t prec and lags 
behind in precession; see Figure JT]). Such a star feels a strong, coherent torque 
from the other stars in the disc, in the opposite direction of its angular momen- 
tum vector, J = \J GM,a(l — e 2 ). As a result, its angular momentum decreases 
in magnitude, which is equivalent to saying that its eccentricity increase fS In 
this way, very high eccentricities can be achieved. Conversely, if a stellar orbit 
is slightly less eccentric than the average, it has a smaller t piec and moves ahead 
in precession, thus experiencing a torque which decreases its eccentricity further 
below the average. 



1 Now, with an even higher eccentricity, the stellar orbit lags even further behind the bulk of the 
disc and the cycle repeats. We emphasize that the precession due to the cusp is a key element 
of this instability. Without it, stars in the eccentric disc would experience torques but in both 
directions and it would not be an unstable configuration. 



Eccentric Disc Instability 
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2.1. Dependency on Background Stellar Cluster 

The stellar orbits at the innermost radii of the disc (i.e, with the smallest a) 
undergo the greatest fractional change in J and hence are most likely to be 
pushed to extreme eccentricities. Consequently, the precession rates of these 
orbits are important. Recalling that t prec oc a 7_3//2 /(e), we perform iV-body 
simulations of an eccentric stellar disc to examine the effect of varying the index 
of the power-law slope of the cusp, 7, on the instability; see Figure ([2]) for 
examples where 7 = 0.5 (shallow cusp) and 7 = 1.75 (steep cusp). The disc has 
a surface density profile S oc r~ 2 , Mdi sc = 10 4 Mq, a semi major axis distribution 
of 0. 05 < a < 0.5 pc and a smooth stellar cusp with M cusp (< lpc) = 0.5 x 1O 6 M 
(see iMadigan et all [2009. for details). For 7 < 3/2, due to the shallowness of 



the cusp, the innermost orbits lag behind in precession and are predisposed to 
being torqued to high e orbits; the instability is not as effective at generating 
low e orbits. Conversely, an increase in 7 (> 3/2) results in the innermost orbits 
precessing faster than the bulk of the stars and hence suppresses the generation 
of the highest e orbits. 



orbits 800 orbits 1600 orbits 









Figure 2. Evolution of the (x and y components of the) eccentricity vectors 
(EVs) of stars within an eccentric disc embedded in a background cusp with 
index 7 = 0.5 (top) and 7 = 1.75 (bottom). Time is in units of the initial 
period of the innermost orbit at a = 0.05pc; final time corresponds to ^0.8 
Myr. The innermost circle shows the initial value of the eccentricities e = 0.6; 
the outermost is e = 1. The colour of the EVs indicate the semi major axes 
of the stars in the disc, going from red (~ 0.05 — 0.1 pc) outwards to 0.5 pc 
(green, cyan and dark blue). The EVs spread out as the stars complete more 
and more orbits, and precess with retrograde motion (clockwise here). Many 
more high eccentricity stars are produced in simulations with a flatter cusp 
(7 = 0.5) as the innermost stars lag behind in precession. In simulations with 
a steep cusp (7 = 1.75), the stars precess more rapidly (t prec ~ 900 orbits) 
and the innermost orbits have lower average eccentricities. 
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2.2. Implications for the Galactic Centre 

A (so far unknown) percentage of stars in a disc should be formed in binary 
systems. These binaries can be propelled to high eccentricities by the ecc entric 
disc instability, pass within the tidal radius of the MBH and be disrupted (jHills! 
119881 : lYu fc Tremainell2003 ). In the GC, this mechanism m ay be responsible for 
producing both hypervel ocity stars (HVS; se e Brown et al.ll2009l. and references 
there i n) and the S-stars (ISchodel et al.ll200"2l ; iGhez et al.l l2005; Eis enhauer et al 



120051 : iGillessen et al.l 120091 ). There are three lines of evidence to support this 
idea. (1) The c usp of late-type stars i n the GC is observed to have a flat surface 
density profile (jBuchholz et al.ll2009l ; [Do et al.ll2009l ; iBartko et al.ll2009l ) which, 
if it traces the underlying mass distribution, is m aximally effectiv e in producing 
high eccentricity orbits with the instability. (2) iLu et al" ( 20091 ) find that the 
observed HVS are consistent with two planar structures, one alined with the 
inner edge of the clockwise disd3- The HVS travel time from the GC to their 
current positions is too large (100 - 200 Myr) for them to have originated in the 
observed you n g disc but an older disc could potentially explain this feature. (3) 
Bartko et al.l ( 20091 ) find a bimodal distribution of eccentricities in the young 



disc. The eccentric disc instability naturally produces a double-peaked eccen- 
tricity profile as the stars are pushed away from the average value. However, 
we must emphasize that a direct comparison to the observations is not possible 
at this time as the most eccentric stars in our simulations are the most highly 
inclined and would not be observed along the plane of the disc. 
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2 The instability preferentially propels stars at the inner edge of the disc to extreme eccentricities. 
These stars are likely to have significant inclinations however which suggests that, in this 
scenario, the HVS will not be tightly confined to a plane. 



